Kinetic Signatures and Intermittent Turbulence in the Solar Wind Plasma 
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A connection between kinetic processes and intermittent turbulence is observed in the solar wind 
plasma using measurements from the Wind spacecraft at 1 AU. In particular, kinetic effects such as 
temperature anisotropy and plasma heating are concentrated near coherent structures, such as cur- 
rent sheets, which are non-uniformly distributed in space. Furthermore, these coherent structures are 
preferentially found in plasma unstable to the mirror and firehose instabilities. The inhomogeneous 
heating in these regions, which is present in both the magnetic field parallel and perpendicular tem- 
perature components, results in protons at least 3-4 times hotter than under typical stable plasma 
conditions. These results offer a new understanding of kinetic processes in a turbulent regime, where 
linear Vlasov theory is not sufficient to explain the inhomogeneous plasma dynamics operating near 
non-Gaussian structures. 



Introduction. — A plasma turbulence cascade [l|, |2| 
might provide the energy needed to heat the lower solar 
corona [3|, accelerate fast and slow wind streams, and ac- 
count for the non-adiabatic expansion of the solar wind 
[j|. However, turbulence at kinetic scales remains an 
unsolved problem and it is not known how fluctuation 
energy is ultimately converted into heat. There are a 
number of kinetic processes, such as linear wave damping 
[5| and pressure-anisotropy instabilities fi], that might 
play a role in this dissipation of the interplanetary cas- 
cade. These are often investigated using linear and quasi- 
linear approximations to the Vlasov-Maxwell equations. 
Other studies have suggested solar wind discontinuities 
are linked to turbulence intermittency, in the sense of co- 
herent structures and non-Gaussian probability distribu- 
tion functions [7| , and are also sites of enhanced temper- 
atures jsl] . Here we ask if kinetic effects are homogeneous 
in space, or concentrated in regions of the turbulent field 
associated with discontinuities. We find a link between 
magnetic coherent structures and kinetic signatures usu- 
ally associated with linear instabilities. 

Kinetic Physics. — The kinetic plasma properties of so- 
lar wind proton populations have been studied in detail 
using in-situ spacecraft measurements |9|. Where sam- 
ples have been subject to many collisional effects [10| . 
the proton velocity distribution functions (VDF) are typ- 
ically isotropic and Maxwellian. However, less collisional 
solar wind plasma usually exhibits anisotropic VDFs with 
respect to the local magnetic field direction, such that 
R = T±/Tu ^ 1. Solar wind temperature anisotropy 
R is studied by investigating the dependence on proton 
parallel beta /3|| = npkBT\\l{B'^ /2^q), which is the ra- 
tio of parallel pressure to total magnetic pressure. These 
parameters cannot assume arbitrary values in the solar 
wind plasma, but are instead confined such that the range 
of observable R is narrowed as /3|| increases. This restric- 
tion could be linked to large deviations from R — 1 which 
can give rise to anisotropy-driven instabilities. These are 
not specific to the solar wind, and can operate in other 



low coUisionality astrophysical plasmas. As the electro- 
magnetic fluctuations associated with these instabilities 
grow, they are able to scatter particles and eventually 
drive the VDF back towards isotropy. From linear Vlasov 
theory, the solar wind plasma can become unstable to the 
mirror and cyclotron instabilities when i? > 1 , and to the 
firehouse instability when R < \ and /3|| > 1. These cal- 
culations were compared with measured distributions of 
R with respect to j3\\ , and the data was found to be best 
constrained by the mirror instability for R > 1 and the 
oblique firehose instability for i? < 1 



mm 



The solar wind plasma is found to migrate with in- 
creasing heliocentric distance towards higher /3|| and to 
values of R approaching unity [13[ . At least part of this 
evolution is better understood by ordering the data ac- 
cording to collisional age r — VppL/vsw, defined as the 
Coulomb proton-proton collision frequency multiplied by 
the transit time from the Sun to 1 AU, which is the num- 
ber of thermalization timescales that elapsed on transit 
from the Sun to the spacecraft [10|. For low collisional 
age plasma, enhancements in power associated with mag- 



netic fluctuations exist near instability thresholds [14 1 
which also correspond to sites of elevated proton heating 
[15| . Indeed, plasma unstable to the mirror or firehose 
instabilities were found to be significantly hotter than 
stable plasma [16J. These associations suggest underly- 
ing physical relationships between heating mechanisms in 
the solar wind and the kinetic physics that emerges from, 
or leads to, the growth of these instabilities. However, 
the details of these relationships are not completely un- 
derstood. This Letter reports results linking instability 
thresholds and elevated temperatures to coherent struc- 
tures dynamically generated by strong MHD turbulence. 

Analysis. — We use around 2.5 x 10^ independent 
plasma and magnetic field measurements from the Wind 
spacecraft, recorded during the interval 1995 Jan. 1 to 
2004 Nov. 20. The Faraday cup instrument of the So- 
lar Wind Experiment (SWE) [17| measures proton den- 
sity Tip, bulk velocity ygw, and proton temperatures. 
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FIG. 1. Plot of median (a) scalar proton temperature Tp, and (b) PVI statistic / over the my ,Tx/r||j -plane. The curves 
indicate theoretical growth rates for the mirror (dashed), cyclotron (dot-dashed), and oblique firehose (dotted) instabilities. 
There is a manifest association between these thresholds, hot plasma, and enhanced PVI. 



These are separated into parallel temperature T|| and 
perpendicular temperature T± by comparison with the 
mean magnetic field from the Magnetic Field Investiga- 
tion (MFI) [18[. Only solar wind data is used, and mea- 
surements either in the magnetosphere or contaminated 
by terrestrial foreshock are removed. We also require the 
uncertainties in the thermal speeds to be < 10% and, to 
avoid Coulomb relaxation effects, the collisional age must 
be T < 0.1. These criteria follow [16|, and only about 28% 
of our original dataset satisfy all these conditions. 

Here we ask if kinetic effects, such as temperature 
anisotropy and heating, are strongly inhomogeneous and 
related to the intermittent character of the turbulent 
magnetic field. A way to find regions of high magnetic 
stress and coherent structures is to identify rapid changes 
in the magnetic field vector: 



AB(i, At) = B(i + At) - Bit) 



(1) 



where B(i) is the magnetic field time series and At is the 
time lag. The fastest available cadence of the plasma data 
defines the lag At = 92 s which, using Taylor's hypoth- 
esis, corresponds to a spatial separation in the plasma 
frame (Ar = —v^wAt) on the scale of inertial range fiuc- 
tuations. We use the partial variance of increments (PVI) 
method to identify coherent (non-Gaussian) structures in 
the solar wind [7|], where the hig hest PVI amplitudes are 
found at the smallest scales [19|: 
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where (• • •) denotes an appropriate time average over 
many correlation times. Events are selected by impos- 
ing thresholds on the I{t) series, leading to a hierarchy 



of coherent structure intensities [8|, |20[ . The PVI statistic 
is constructed such that (/^) = 1, (I^) is related to the 
kurtosis of the magnetic field increments, and powers of 
(/) scale in a manner connected with familiar diagnostics 
of intermittency. 

Results. — The scalar proton temperature Tp and PVI 
statistic were each divided into 50 x 50 grids of logarith- 
mically spaced bins in the (/?|| , R) plane, and the median 
values in each bin containing at least 50 measurements 
were computed. Figure la is (by construction) similar to 
Fig. 2 in |16j . and solar wind observations are confined 
in the {(3i\, R) plane in a manner consistent with previ- 
ous studies [e.g. llJ-llJ|. The highest Tp values occur in 
regions typically associated with linear instability thresh- 
olds, and even at high /3|| these areas are around 3-4 times 
hotter than the nearby R = I plasma. However, the 
instabilities themselves are not thought to significantly 
heat the plasma [21|, and so the heating mechanism has 
not been conclusively identified. 

Insight into the heating mechanisms contributing to 
these observations can be obtained by considering the 
PVI statistic. Figure lb shows the median PVI values in 
the (/3||, R) plane. The elevated / values, which corre- 
spond to an increased likelihood of finding coherent mag- 
netic structures such as current sheets within the observa- 
tions, occur in almost exactly the same regions associated 
with enhancements in Tp. This result is completely con- 
sistent with earlier findings [20] that samples of stronger 
PVI events systematically produce conditional probabil- 
ity distributions with higher mean proton temperature. 
Indeed, the same procedure also produces conditional dis- 
tributions of electron temperature and electron heat fiux, 
all with increased average values. However, it was not 



previously known that the most intense PVI events are 
preferentially found in plasma with high values of tem- 
perature anisotropy. This suggests coherent structures 
dynamically generated by MHD turbulence could be re- 
sponsible for the anisotropic solar wind heating shown in 
Fig. la and observed by |15llla|. 

The relationship between PVI and Tp is further inves- 
tigated by dividing both into a 25 x 25 grid of logarith- 
mically spaced bins, and computing the joint probability 
distribution p{I, Tp) by dividing the number of observa- 
tions in each bin n by the total number of observations N. 
Figure 2a shows the proton temperature as a function of 
PVI, and the white dots are the most likely value of Tp in 
each PVI bin. There is a clear dependence of proton tem- 
perature on PVI, and the most significant enhancements 
in Tp are associated with PVI events above a threshold of 
around / > 3. These events correspond to non-Gaussian 
coherent current sheet structures [7[ , where the most in- 
tense are candidate magnetic reconnection sites [l9|. 
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FIG. 2. Joint probability distribution between (a) proton 
temperature and PVI, where the most likely values of Tp 
in each PVI bin are represented by white dots. There is a 
clear link between enhanced Tp and strong PVI events. The 
most likely (b) temperature anisotropy value for each PVI 
bin, where Tl (Ts) is the largest (smallest) of T± and Ty. 
The most intense PVI events are associated with the greatest 
departures from temperature isotropy. 

A link between the most intense PVI events and en- 
hancements in Tp has been established [8|, |20[. However, 
it is unknown if a similar link exists between PVI and 
temperature anisotropy in the solar wind. Here we de- 
fine Tl and Ts such that Tl/Ts = R when R > 1 
and Tl/Ts = R'^ when R < I. Figure 2b plots the 
most likely values of Tl/Ts, and the associated stan- 
dard errors, for each of the 25 logarithmically spaced PVI 
bins. The temperature anisotropy remains almost con- 
stant until around 1 = 3 when it starts to increase with 
stronger PVI events. This is consistent with tempera- 
ture anisotropy being inhomogeneous and concentrated 
in the vicinity of coherent structures. These results are 



corroborated by Vlasov-Maxwell simulations, where tem- 
perature anisotropy is found to be enhanced in regions of 
strong magnetic stress between magnetic vortices 22 1 . 

Temperature anisotropy near coherent structures could 
manifest through the preferential heating or cooling of 
protons either parallel or perpendicular to the mean mag- 
netic field. Figure 3 shows PDFs of T|| and T± condi- 
tioned on PVI value, where / > is the entire dataset, 
/ > 1 removes low value fluctuations, / > 3 only retains 
non-Gaussian structures, and / > 5 contains only the 
most highly inhomogeneous structures including current 
sheets. Plasma with the greatest proportion of intense co- 
herent structures has the highest most probable parallel 
and perpendicular temperatures. However, this heating 
appears in anisotropic proportions in the vicinity of the 
most non-Gaussian structures. Therefore, this suggest 
the presence of dynamical heating mechanisms within or 
nearby these coherent structures that can heat and ac- 
celerate protons both parallel and perpendicular to the 
magnetic field direction. 

Disscussion. — The results detailed here suggest dis- 
tinctive kinetic signatures observed in solar wind plasma 
are associated statistically with coherent magnetic struc- 
tures such as current sheets, which are connected to 
the intermittency properties of MHD turbulence. These 
kinetic effects include plasma heating and temperature 
anisotropy. The elevated proton temperature is due to 
broadening of the underlying VDF, which could be either 
a consequence of energy dissipation or particle energiza- 
tion, or a combination of both. Furthermore, plasma at 
the extremes of the accessible (/3|| , R) parameter space 
which is bounded by linear instability thresholds is not 
only likely to be hotter, but also likely to contain stronger 
magnetic discontinuities (as measured by higher PVI). 
Indeed, the pattern of median PVI values in this plane 
closely resembles the contours of median temperature. 

We are led to the conclusion that kinetic effects in 
the solar wind plasma are inhomogeneous, being con- 
centrated at the most intense coherent structures which 
are related to the intermittent properties of turbulence. 
However, the exact nature of this relationship between 
plasma turbulence and kinetic effects is unclear. It is 
known that MHD turbulence dynamically generates co- 
herent structures, and these are sites of enhanced heating 
as in Fig. 2a Isl. 12011 and elevated temperature anisotropy 
as in Fig. 2b [22|. This turbulence produced tempera- 
ture anisotropy could then move the plasma to marginal 
instability, driving the growth of the firehose and mir- 
ror instabilities at the appropriate boundaries. Alter- 
natively, the growth of these instabilities might result 
in the generation of coherent structures which heat the 
surrounding plasma anisotropically, thus driving it back 
to marginal stability. Indeed, there could be a complex 
feedback mechanism underlying the interaction between 
turbulence and these kinetic phenomena. However, the 
limiting regions of the {/3h , R) plane seem to be roughly 
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FIG. 3. PDFs of the (a) parallel and (b) perpendicular proton temperatures, where each corresponds to a different range of 
PVI values. In both cases, the strongest PVI events are associated with elevations in temperature. 



as well correlated to enhancements in the PVI statistic as 
they are to hnear Vlasov instability thresholds. There- 
fore, consideration must be given to the possibility that 
an explanation to our observations need not include ef- 
fects deriving from linear Vlasov theory. 

There is evidence supporting the link between turbu- 
lence generated coherent structures and kinetic phenom- 
ena such as temperature anisotropy and plasma heating. 
This motivates further study to understand the detailed 
plasma physics underlying these connections. In partic- 
ular, examination of dynamical activity associated with 
current sheets such as magnetic reconnection. For exam- 
ple, test particle scattering in MHD simulations shows 
most energetic protons are accelerated in or near cur- 
rent channels due to interaction with the nearby inho- 
mogeneous electric field. These particles are not only 
accelerated to suprathermal speeds, but rapidly form 
anisotropic velocity distributions with speeds perpendic- 
ular to the mean magnetic field enhanced relative to 
speeds parallel [23|. There have also been indications 
of anisotropic temperatures near current sheets in 2.5D 
kinetic hybrid [24] and hybrid Vlasov 22| plasma tur- 
bulence simulations. Hence, there is strong impetus to 
understand dissipative and heating processes that oper- 
ate non-uniformly in space. This should supplement the 
standard approach that focuses on linear models of ki- 
netic dynamics within a uniform homogeneous plasma. 
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